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ABSTRACT 

in 

o. 

^D . We have investigated the young stellar population in and around SFO 38, one of 

the massive globules located in the northern part of the galactic H ii region IC 1396, 
using the Spitzer IRAC and MIPS observations (3.6 to 24 fira) and followed up with 

^ ' ground based optical photometric and spectroscopic observations. Based on the IRAC 

d • and MIPS colors and Ha emission we identify ~ 45 Young Stellar Objects (Classes 

O/I/II) and 13 probable Pre Main Sequence candidates. We derive the spectral types 
(mostly K- and M-type stars), effective temperatures and individual extinction of the 
relatively bright and optically visible Class II objects. Most of the Class II objects 
show variable Ha emission as well as optical and near-infrared photometric variability 
which confirm their "youth" . Based on optical photometry and theoretical isochrones, 
we estimate the spread in stellar ages to be between 1-8 Myr with a median age of 
3 Myr and a mass distribution of 0.3-2.2 Mq with a median value around 0.5 Mq. 
Using the width of the Ha emission line measured at 10% peak intensity, we derive the 
mass accretion rates of individual objects to be between 10^^*^ to 10^^ Mq/jt. From 
the continuum-subtracted Ha line image, we find that the Ha emission of the globule 



is not spatially symmetric with respect to the O type ionizing star HD 206267 and 
the interstellar extinction towards the globule is also anomalous. We clearly detect an 
enhanced concentration of YSOs closer to the southern rim of SFO 38 and identify an 
evolutionary sequence of YSOs from the rim to the dense core of the cloud, with most 
of the Class II objects located at the bright rim. The YSOs appear to be aligned along 
two different directions towards the 06.5V type star HD 206267 and the BOV type star 
HD 206773. This is consistent with the Radiation Driven Implosion (RDI) model for 
triggered star formation. Further the apparent speed of sequential star formation is 
consistent with the speed of propagation of shocks in dense globules as derived from 
numerical simulations of RDI. 

Subject headings: stars: low mass, pre-main-sequence, ISM: clouds and H ii regions: 
individual: SFO 38, IC 1396 



1. Introduction 



Observational and theoretical studies of star formation over the last decade have increasingly 
strengthened the idea that massive young stars play an important role in triggering the formation 
of subsequent generations of stars. The triggers of star formation which are typically in the form 
of wind and radiation from the massive stars or expansion of H ii regions essentially involve the 
induction of compression of a molecular cloud externally. Among several proposed theories, two 
models have gathered sufficient observational support in order to be regarded as the most plau- 
sible models for triggered star formation. In the firs t model known as the Collect and Collapse 
(jElmegreen &: Ladalll977l : iHosokawa &: Inutsukall2006l ) model, an expanding H ii region sweeps up 
material into a dense bordering layer between the H ii region and the molecular cloud. The com- 
pressed shell of dust and gas undergoes fragmentation and gravitational collapse to form new stars. 
Observational evidence supporting the "collect and collapse" model includ es detection of a dense , 



fragmented shell of gas with new ly formed stars surrounding the H ii region (JDeharveng et al. 



2009 



2003 



Zavagno et al.ll2006l . 120071 ). The collect and collapse model is however primarily to explain the 



triggered formation of massive stars. An equally likely model based on Radiation Driven Implosion 
(RDI), involves the creation of shock front at the surface of molecular clouds due to photoevapo- 
ration from the surfaces of molecular clouds exposed to the UV radiation of nearby massive young 
stellar clusters. The enhanced inward pressure triggers the formation of new protostellar cores or 
compresses pre-existing ones to form a new generation of stars. Observational evidence for the 
RDI model includes typical sp atial distribution and gradient in evolutionary phases of young stel- 
lar objects (jlkeda et al.ll2008l ). There is however considerable ambiguity regarding whether RDI 
actually induces the gravitational collapse or merely exposes the young stars by photoevaporation, 
since these two scenarios are observationally indistinguishable. Thus additional tests to ascertain 
the temporal feasibility of shock waves from nearby OB associations actually triggering sequential 
star formation are required to strengthen the arguments in favor of RDI. 



Bright-rimmed clouds (BRCs) are isolated molecular clouds located on the edges of evolved 
H II regions, which are found to have many signposts of sta r forn iation, i.e., IRAS poin t sources, 



molecular outflows, HH objects etc. ISugitani et al.l (jl999l ) and ISugitani et al.l (J2000|) observed 



a sample of 89 BRCs and found an excess in the luminosities and luminosity-to-cloud ratios of 
embedded IRAS sources when compared with sources in isolated globules, which is an indication 
of enhanced star formation in the BRCs. Young stellar objects seemingly aligned along the axis 
towards the ionizing cluster were detected via near- infrared (NIR) imaging of 44 of these BRCs. 
There is also evidence for an age gradient, with older stellar objects closer to the OB cluster and 
the younger objects well inside the globules, aligned with the IRAS sources. These results are 
consistent with sequential formation of stars while the shock front advances further and further 
into the molecular cloud. Thus BRCs are ideal objects for studying and verifying the different 
models of triggered star formation. 



The H II region IC 1396, powered by the 06.5V type (|Sticklandlll995l l star HP 206267 i n the 



Trumpler 37 cluster, appears to sweep up a molecular ring of radius 12 pc (IPatel et al. 



is surrounded by 11 Brigh t-Rimmed clouds with embedded IRAS point sources (ISugitani fc Ogura 



1994 



1998) and 



Sugitani et al.lll99ll ). SFO 38, also known as IC 1396N, is located at a projected distance of 



~11 pc to the north of HD 206267 with the bright rim corresponding to the ionization f ront facing 
the 6.5V star. SFO 38 is located in the Cep 0B2 association at a distance of 750 pc (JMatthews 
19791 ). More recent Hip parcos parallax based measurements estimate the distance to Cep 0B2 
association to be 615 pc (jde Zeeuw et al.lll999l ). However in order to make our results comparable 
to most published results, in this paper we adopt a distance of 750 pc for SFO 38. The derived 
luminosities may therefore be overestimated by more than 40% and physical sizes by 20%. The 
region is associat ed with IRAS 21391+ 5802, a very young intermediate-mass objec t, with a lumi- 



nosit y of 235 Lq (ISaraceno et al.lll996l ). which powers an exte nded bipolar outflow (jSugitani et al 



19891 ). Based on millimeter observations iBeltran et al.l ([2002) resolved IRAS 21391+5802 into an 
intermediate-mass sourc e named BIMA 2 surro unded by two less massive and smaller objects, 
BIMA 1 and BIMA 3. IValdettaro et al.l (|2005l ) detected II2O maseremission at 2 2 GHz to- 
wards SFO 38 which is consistent with an intermediate-mass object. iNeri et al.l (|2007l ) used still 
higher angular resolution millimeter interferometric observations to reveal that the intermediate- 
mass protostar BIMA 2 itself consists of multiple compa ct sources. The gas emission surround 



ing I RAS 21391+5802 traces different molecular outflows (jCodella et al. 



2001 



Beltran et al. 



2002, 



20041 ). NIR images of the region have also revealed a number of small scale molecular hydro- 



gen and Herbig-Haro (HH) f lows (INisini et al. 



Caratti o Garatti et al. 



2006 



Beltran et al. 



2001 



200 



i" 



Sugitani et al.l I2OO2I : iReipurth et al. 



2003 



These observations confirm the on going star 



formation in the dense core of SFO 38. 



Getman et al.l (|2007l ) identified the mid-infrared (MIR) (Spitzer) counterparts of X-ray sources 



detected in Chandra observations and found that the young stellar objects (YSOs) are oriented 
toward the ionizing star HD 206267 and show an age gr adient consistent with the RDI model for 
triggered star formation. Recently iBeltran et al.l (J2009l ) studied the YSO population of SFO 38 



by obtaining deep J, H, K^ broadband images and deep high-angular resolution observations in 
the H2 narrowband filter. iBeltran et alj (|2009l ) did not detect any clear NIR exce ss close to the 



rim, which is in contrast to the age 1 Myr of the stellar population concluded by iGetman et al 



(|2007|). iBeltran et all ((20091) suggest that the YSOs closer to HD 206267 could actually be younger 
than the Class II objects they appear to be since their circumstellar environment w ere disrupted 
completely by intense UV radiation field from the OB association. IBeltran et al.l thus suggest 
that in general the apparent age sequence seen close to OB associations need not necessarily be 
the actual evolutionary sequence and in the case of SFO 38 there is no concrete evidence of star 
formation by the nearby OB association. 

In this paper we address the controversial issue of the distribution of YSOs by making use of 
the extremely sensitive Spitzer IRAC (3.6 to 8.0 iJ,m) k. MIPS (24 /xm) dat a in order to p r obe a ll 
embedded YSOs and protostars in SFO 38. This thus extends the work by iGetman et al.l ()2007l ). 
who had considered the 3.6-5.8 //m characteristics only of those MIR sour ces from which X-ray 
emission had been detected. Our approach also complements the work by IBeltran et al.l (120091 ) 
which makes use of only the NIR observations to decipher any age gradient of the YSOs. In this 
work we shall show that most of the Class O/I sources identified in the MIR are either barely 
detected (with large uncertainties in fluxes) or not detected at all in the NIR and in the X-rays. 
In this paper we have completed the multiwavelength overview of SFO 38 by obtaining optical 
broadband images in B VRI bands and narrow band Ha image to understand the detailed structure 
of the bright-rim. Further in order to confirm the ^^youth" of the YSOs identified in the MIR we 
have made use of medium resolution spectroscopy of the optically visible YSOs in Ha line emission. 
Observation details of Spitzer and optical 5 Vi?/ imaging and medium resolution spectroscopy are 
presented in § [21 We describe the selection of YSO candidates based on MIR color-color diagram 
(CCD) and the detailed analysis of optical-MIR data in section §|3l In §|l|we discuss the probable 
star formation scenarios and in § [5] we summarize the results. 



2. Observations and Data Reduction 

2.1. Spitzer IRAC & MIPS observations 

We have extracted IRAC (3.6, 4.5, 5.8, & 8 ^m) and MIPS (24 & 70 fim) observations from the 
Spitzer Space Observatory archive (Program ID 30050: Star Formation in Bright Rimmed Clouds 
by Fazio et al.). The IRAC data were taken in the High Dynamic Range (HDR) mode using a single 
AOR (Astronomical Observation Request) with a five-point dither pattern. We have processed both 
the short (0.6 sec) and the long (12 sec) integration Basic Calibrated Data (BCD) frames in each 
channel using the Artifact mitigation software developed by Sean Carey lj and created mosaics 
using MOPEX. In all IRAC bands we detect point-sources down to 70 /iJy. We have also created 



^ |http : //web. ipac.caltech.edu/steLff/caj:ey/irac_aj:tif acts/index .html I 



mosaics of the MIPS 24 and 70 /iin BCDs using MOPEX. Both data sets are of reasonably good 
quaUty. The 70 /im-image, due to the lower angular resolution and smaller area mapped, shows 
only one bright point-like source embedded in the globule. 

We have carried out multiframe PSF photometry using the tool APEX developed by Spitzer 
Science Center(SSC) on all the Spitzer IRAC images and on the MIPS images. The long integration 
3.6 and 4.5 ^m IRAC images show signs of saturation on bright stars. For the saturated sources, 
the photometry derived from the short integration images were used. It is difficult to disentangle 
the sources in regions of strong emission from the associated Photon Dominated Regions (PDRs) 
from the surrounding clouds and to derive their accurate photometry. We have used a combination 
of automated routines and eye-inspection to detect sources and extract photometry of these sources 
from the IRAC and MIPS images. For sources, which APEX failed to detect automatically at one 
or several wavelengths, we have used the user list option in APEX to supply the coordinates for 
the source to successfully derive a PSF fit. This enabled us to derive photometry for every source 
which we could visually identify on any image. Following the IRAC and MIPS Data Handbooks 
we have adopted the zero-points for conversion between flux densities and magnitudes to be 280.9, 
179.7, 115.0, 64.1 and 7.14 Jys in the 3.6, 4.5, 5.8, 8.0, and 24 /um bands, respectively. 



2.2. Optical BFi?/ Photometry 

SFO 38 was observed on 17th July, 23rd August and 15th September 2009 using Bessell's broad 
band filters VR (300x8 s, 300x8 s), 57/(900x8 s, 300x8 s, 150x8 s) and 5/(900x8 s, 120x12 s) 
respectively and on 4th November 2008 using the narrow band Ha (900 s) with the Himalayan Faint 
Object Spectrograph Camera (HFOSC) mounted on the 2m Himalayan Chandra Telescope (HCT) 
of the Indian Astronomical Observatory, Hanle, Indiao- HFOSC uses the central 2Kx2K region of 
the 2Kx4K CCD in imaging mode and covers a field of 10 xlO with a plate sca l e of (y .'296 pixel~^ . 



The nights were photometric with an average seeing of 1'.'5 to 1'.'8. iLandolti (119921 ) photometric 
standard stars were observed on all the nights to calibrate the target stars. Data were reduced using 
standard tasks available within Image Reduction and Analysis Facility (IRAFJfj. Bias subtracted, 
flat field corrected and aligned frames were combined to make the master frame for each filter. 
Astrometric calibration was applied to the master frames using the IDL procedure STARASTT of 
IDL Astronomy User's Librar}Q. Photometric magnitudes were calculated by aperture photometry 
with the optimal aperture adopted as the radius where the difference in magnitudes between two 
consecutive apertures is less than 1%. The limiting magnitude is defined as the magnitude at which 
the mean magnitude error of the star becomes 0.1 mag which implies a 10a detection corresponds 



•^ http : //www. iiap . res . in/centers/iao | 

^The IRAF software is distributed by the National Optical Astronomy Observatory under contact with National 
Science Foundation. |http : //iraf . noao . edu/] 

^ [http : //idlastro . gsf c . nasa . gov/] 



to a signal-to-noise ratio (S/N) of 10. The limiting magnitudes were V: 21, R: 21 on 17th July , 
B: 22, V: 22, I: 19.5 on 23rd August and B: 22.2, I: 20.2 on 15th September. Aperture corrections 
for each frame were derived from the bright and isolated stars and applied to the faint stars. The 
standard deviation of residuals of observed and transformed magnitudes and colors of the standard 
stars are within the range of 0.01-0.02 mag. We used the broadband R filter image for con t inuum 



subtraction from the narrow band Ha image of SFO 38. Based on the recipe of IWalleii (|l99Cl ). 
the point spread functions of the R-band and the narrow band Ha images were matched and the 
images were scaled before subtraction to get the Ha emission line image of SFO 38. 



2.3. Medium Resolution Spectroscopy 

Medium resolution (~7 A) spectra were obtained for relatively bright and optically visible 
YSOs during July-November of both 2008 and 2009 in the wavelength range (5200-9200 A) with 
the Himalayan Faint Object Spectrograph Camera (HFOSC) mounted on the Himalayan Chandra 
Telescope (HCT). Typical exposure time of each spectrum was 3600 s. The spectra were bias 
subtracted, flat fielded and the one-dimensional spectra were extracted using the standard tasks 
of Image Reduction and Analysis Facility (IRAF). The arc lamp spectra of FeNe were used for 
wavelength calibration. We further used the strong night sky emission lines e.g. [01] A5577 A, 
A6300 A lines and rotational and vibrational bands of OH in the red region of the individual target 
spectrum to improve the wavelength calibration and achieved an accuracy of it 0.5 A for each 
target star. 



3. Results and Analysis 

3.1. Mid-Infrared (3.6 to 24 /xm) view of SFO 38 

Figured] presents the three-color image of SFO 38 using the IRAC 3.6, 8.0 and the MIPS 24 /im 
bands. The image shows that the emission from the front side of the globule is dominated by strong 
polyaromatic hydrocarbon emission (PAH-emission) , and rather clumpy. This could in part be due 
to Rayleigh- Taylor instabilities at the ionization front, but it is also due to the outflows from the 
young stars. The cometary shape of the globule is quite apparent with the tip pointed towards the 
O star HD 206267 and the eastern edge appears to be more abruptly truncated to the south. This 
indicates that the south-eastern edge of the globule has most likely experienced intense ionizing 
radiation field, which has eroded the cloud material more than on the western side. Comparison 
with the Ha emission line image (Fig. [2]) also confirms the presence of a more pronounced ionization 
front to the south-east and a rather tenuous layer of ionized gas to the west of SFO 38. 

The aim of this paper is to identify the young and embedded stellar population in SFO 38. 
While the 3.6 and 4.5 ^m IRAC data is particularly sensitive, a large number of stellar objects in 



addition to the YSOs also detected in these bands. We have thus used all sources detected either 
at 5.8 or at 8.0 fim or at both wavelengths to create a list of MIR sources in SFO 38. Photometry 
for all these sources, if detected, are then extracted from the 3.6 &: 4.5 //m images. We thus detect 
98, 106, 106 and 98 sources respectively in the IRAC 3.6, 4.5, 5.8 and 8.0 //m images and 14 and 
1 sources respectively in the MIPS 24 &: 70 fiui images. All the sources detected at 24 /im are 
found to have been detected in the IRAC bands. In all we identify a total of 110 MIR sources. 
We have crosscorrelated the MIR sources with the NIR sources in our mapp ed region from the 
2M ASS point source cat alogue, as well as the sources detected in the NIR by iNisini et al.l (J200ll ) 
and iBeltran et al.l (|2009l ). We have used the following association radii : 1" for the IRAC images, 
2'.'5 for the MIPS 24 /im image and 2" for NIR data. Table [J gives the coordinates of the 110 sources 
identified in SFO 38 together with their NIR magnitudes, Spitzer IRAC and MIPS flux densities 
and a preliminary classification based on selected color-color plots and other criteria described in 
Sec. 13.21 In Table [1] we have used the prefix SF038, but throughout the paper we refer to them 
as MIR-nn, where nn is the number of the source. Out of the 110 sources 80 were found to have 
NIR counterparts. However 10 sources o ut of the 80 sources with NIR counterparts were found to 
coincide with the H2 knots identified by lBeltran et al.l (|2009l ) . 



Additionally in order to identify the stellar and pre main sequence (PMS) stars we used pho- 
tometry extracted from the IRAC 3.6 and 4.5 /um long integration images. The 3.6 and 4.5 iim 
images are more sensitive, have a cleaner PSF, and appear to be less affected by nebular emission 
than the 8 //m image. We identified 161 additional sources, which are detected in both the 3.6 and 
4.5 /um wavebands. Of these sources 113 are found to have NIR counterparts, and 6 sources are 
coincident with H2 emission knots. 



3.2. Classification of YSOs based on Near- and Mid-infrared colors 



In regions of high obscuration like SFO 38, which has a visual extinction of up to 20 mag, 
mid-infrared color indices provid e an ideal tool for id entifying YSOs and classify them according 
to their phases of PMS evolution (JGetman et al.ll2007l ). Here we use several color indices generated 
using the Spitzer IRAC and MIPS data to characterize the evolutionary stages of the MIR sources 
detected in SFO 38. 



Stern et al.l (120051 ) demonstrated that (a) normal star-forming galaxies and narrow-line AGNs 



with increasing 5.8 and 8.0 /^m and (b) broad-line AGNs with red, nonstellar SEDs, result in colors 
which are very similar to bona-fide YSOs. Thus prior to identifying PMS sources based on the 
color in dices it is necessary to ins pect deep IRAC images for contamination due to the extragalactic 



sources. 



Gutermuth et al.l (|2009l ) have extensively discussed the criteria for identifying such extra- 



galactic objects in the IRAC color-color diagrams. Using the criteria given in the appendix of their 
paper, we identify MIR-2, 17, 68, 70, 75, 81, 84, 89, 101 and 103 as AGNs and MIR-77, 81, 93 and 
103 as galaxies strongly emitting in the PAH bands. Of these sources, MIR-17, 89 and 103 satisfy 
most of the criteria for extragalactic objects. The profiles of sources in the i?-band is typically 



quite a good indicator of whether the source is actually extragalactic or stellar. However all the 
sources in SFO 38 which satisfy the criteria mentioned above are extremely faint so that with the 
exception of MIR-2, these sources though visible are not detected clearly enough in the -R-band 
so that the intensity profile can be studied. Further, none of the sources is detected in the MIPS 
24 fiui image either. The -R-band profile of MIR-2 appears to be stellar with a faint companion 
on one side and MIR-84 is associated with an Ha emission star, so we do not identify both these 
sources as extragalactic. In the absence of any conclusive evidence we presently assume the sources 
identified by the empirical criteria as above to be extragalactic. 

The 4.5 ^m IRAC band of Spitzer covers many of the H2 pure rotational lines, as a result of 
which this band records a large excess of emission at the position of the unresolved blobs of shocked 
emission from high velocity outflows which i nteract with the molecular cloud. Using the empirical 
criteria outlined by iGutermuth et al.l (120091 ) here we flag the sources MIR-24, 27, 53, 65 and 98 
as possible fe a tures due to shocked H2 emission. Based on 2.12 //m H2 integrated line emission 



Beltran et al.l (120091 ) identifled as many as 97 emission features possibly created by H2 excitation 



by shocks driven by outflows powered by YSOs. We find that 10 out of the 110 MIR sources that 
we have identifled are associated with such H2 knots. These include the sources MIR-20, 24, 27, 
35, 53, 65, 74, 75, 98 and 105. Thus the color indices for 5 out of the 10 sources are consistent 
with the observed association with H2 knots. The source MIR- 75 which also satisfied the criteria 
for being an extragalactic object is thus identified as a H2 knot. 

Figure [3] presents color-color diagrams of NIR, IRAC and MIPS 24 /xm sources detected in the 
SFO 38 globule. We have used several criteria (shown as dashed lines and boxes in Figure [3]) to 
identify the potential YSOs using these color-color diagrams. 

The most stringent classification scheme uses the IRAC ([3.6]-[5.8]) colors and the [8]-[24] 
IRAC and MIPS color. At 24 fiui the reddening d ue to extinction is sin all and the photospheric 
colors are very close to zero for all spectral types (JMuzerolle et al.ll2004l ). Therefore the [8]-[24] 
color is very sensitive to infrared excess, but of course not all young stars are bright enough to be 
detected at 24 fan. Using this color-color diagram we find fourteen sources with infrared excess 
(Figure El (Ze/i)). We identify MIR-34, 36, 48, 50, 54, 55 and 59 as Class O/I sources, while MIR- 
11, 31, 38, 45, 82 and 86 are within the Class II regime and MIR-29 which has a large [8.0]-[24] 
color excess and very small value for the [3.6]- [5.8] color is identified as a transition disk object. 



J0rgensen et al.1 (|2007l ) predicted sources with [3.6]-[4.5]>1.0 and [8]-[24] > 4.5 to be YSOs driving 



molecular outfiows. We find that the sources MIR-34, 48, 50, 54, and 59, all of which have been 
identified as Class O/I objects also satisfy the criteria for being sources driving outflows. 

Figuie^ (middle) presents the [3.6]-[4.5] vs [5.8]-[8.0] color-color plot for the sources detected 
in all the IRAC bands. Sources with the colors of stellar photospheres are centered at ([3.6]- 
[4.5],[5.8]-[8.0])=(0,0) and include foreground and background stars as well as diskless (Class III) 
pre-main sequence stars. Sources which satisfy the criteri a — 0.1 < [3 . 6] -[4. 5] < 0.1 and 0.2 < [5.8]- 



[8.0] < 1.1 are classified as the transitional disk objects (|Fang et al.l 12009. ) which also do not have 



any NIR color excess. The box outlined in Figure [3] (middle), defines the location of Class II objects 
(JMegeath et al.ll2004l : lAllen et al.ll2004l'). i.e. source s whos e colors can be explained by young, low- 
mass stars surrounded by disks. iHartmann et al.l (|2005l ) have shown from their observations of 
young stars in the Taurus- Auriga complex that Class O/I protostars require [3.6]-[4.5]> 0.7 and 
[5.8]-[8.0]> 0.7. With the exception of MIR-29, the classes derived for the sources identified in the 
MIPS 24 /um band and plotted in Fig. El^/e/i) are completely consistent with the classes derived from 
the [3.6]-[4.5] vs [5.8]-[8.0] color-color plot. We identify 10 of the MIR sources as being associated 
with H2 knots, this includes MIR-75, which otherwise satisfies the criteria for extragalactic objects. 
We also plot NIR and Spitzer MIR color-color diagram of the identified YSOs and transitional Disk 
objects in Figure [3] (ri(7/i^). 

Table [T] presents a summary of classification of the MIR sources based on the color-color 
diagrams presented here. We have additionally looked into the [3.6]-[5.8] vs [4.5]-[8.0] color-color 
diagram and found that the Classes derived from this diagram are consistent with the classes derived 
based on the color-color diagrams presented here. For most sources it was possible to arrive at a 
unique Class based on all three color-color diagrams. For all sources the Class derived based on at 
least two color-color diagrams match. 

Use of the more sensitive long integration 3.6 and 4.5 //m images in combination with 2MASS 
data (following a classification scheme involving the Ks-[3.()\ and [3.6]-[4.5] colors) yields an addi- 
tional 113 sources, 11 of which are PMS objects. Since the long integration 3.6 and 4.5 fim images 
go deeper and are less affected by extinction than the 2MASS survey, we also checked how many 
of these IRAC sources which have no 2MASS associations are likely to be PMS s tars. Among the 



48 rem aining sources which have no NIR counterparts, following the arguments of lMookeriea et al. 



(J2OO9I ) if we use the color criterion, [3.6]- [4.5] > 0.2, we obtain 8 s ources. However 6 o f these 8 
sources are actually associated with H2 emission knots as detected bv lBeltran et al.l (|2009l ). Table [2] 
lists the 13 PMS candidates identified on the basis of NIR and 3.6 and 4.5 fim magnitudes. 

Excluding the H2 knots and the possible extragalactic objects we identify 10 Class O/I sources, 
3 Class I sources, 13 Class I/II sources (occupying the top- left and bottom-right corners), 14 Class 
II sources and 19 transitional disk objects in SFO 38. Thus we find a total of 40 YSOs (Classes 
O/I/II) and 13 PMS candidates based on MIR color indices. We note that 13 Class I/II objects 
and 19 Transitional Disk objects are ~ 30% of the selected YSOs and show noticeable color excess 
only beyond 6 /um. 



3.3. Ha emission of SFO 38 and spatial distribution of YSOs 



Sugitani et al.l (|l99lh classified SFO 38 as a 5 type BRC powered by 06.5 star HD 206267. 



The bright rim structure is quite prominent in the continuum subtracted Ha line image of SFO 38 
in Figure [2j The Ha emission nebulosity of the bright-rim of the globule is asymmetric with 
respect to the geometric axis of the globule as well as the direction of the adopted ionizing source 



10 



HD 206267. The ionized rim appears to be brighter towards the east and rather diffuse towards the 
western edge, facing HD 206267. The Ha intensity distribution of the diffuse material in between 
the bright-rim and the ionizing star also show similar kind of asymmetric distribution as that of 
the bright-rim. We, therefore, searched for other massive O and B type stars in the surrounding 
region that can p roduce the relatiyely b righter Ha emission on the left side of the globule. We 
find a BOV type (jJaschek fc Egret! Il982l ) star HD 206773 which is situated on the symmetry axis 
of the left side of the bright-rim at a projected distance of ~ 7.6 pc from SFO 38. Ide Zeeuw et al. 



(j 19991 ) tagged HD 206773 as a member of Cep 0B2 association. The distance towar ds the star is 



betw een 550-1000 pc, adopting the parallax measurement and the associated error (jvan Leeuwen 
20071 ) which is consistent with the dis tance of SFO 38. Th e number of ionizing photons emitted 
per second by a BOV star is ~ XO"^^"^ ( Sternberg et al.l 120031 ) and using n// ~ 10^ m~^ for number 
density of hydrogen in the interior of the H II region IC 1396, the Stromgren radius would be ~ 
8.5 pc which suggests that HD 206773 could be a potential ionizing source for SFO 38. Fig. [2]shows 
the distribution of the YSOs and the directions toward HD 206267 and HD 206773 on the Ha line 
image. We find that most of Ha emitting Class II objects lie within the bright-rim. Class I/II 
objects are situated behind the Class II objects in the intermediate region between the bright-rim 
and the dense portion of the globule. All the Class O/I and Class I objects are situated well within 
the dense core part of the globule. There is no Class O/I to Class I/II YSO which is situated before 
the ionized rim. The YSOs appear to have elongated distribution to the right and left sides of the 
bright rim with two different axes of elongation. 



3.4. Additional YSO candidates based on Ha emission 



Ha emission surveys are very useful to identify low mass (< 2 M©) Classical T Tauri stars 
(CTTS) in young star forming regions as the Ha emission is due to the accretion of circumstellar 
material in YSOs. These observations also detect a fraction of Weak-lined T Tauri stars with weak 
Ha emission which is believed to arise from enhanced chromospheric activities. Based on Ha emis- 
sion, we searched for the additional YSO candidates available in the literature outside the field of 
view of Spitzer observations to complement the list of YSOs around SFO 38. Many of the candi- 
date YSOs in and around SFO 38 were found during Ha objective prism surv ey by Ogu ra et alj 
(l2002l'):llkeda et all pOod ). A total of 14 out of the 17 Ha emission line sources (|Ogura et al..,2002i : 



Ikeda et al. 



20081 ) are situated within the field of view of IRAC. We consider the rest of the Ha 



emission line sources viz., [OSP2002] BRC 38 1, [OSP2002] BRC 38 16, and [OSP2002] BRC 38 17 
as YSO candidates. We also find 4 Ha emission line sources near SFO 38 in the catalogue of point- 
source Ha emission-line object s selected froin the INT/WFC Photometric Ha Survey (IPHAS) 
of the northern Galactic plane (jwitham et al.ll2008l ). We identify IPHAS J214027.32+581421.3 as 
[OSP2002] BRC 38 2, IPHAS J214036.90+581437.9 as [OSP2002] BRC 38 6, IPHAS J214041.23-h581158.5 
as [GFG2007] 64. Thus, we also consider [GFG2007] 64 and IPHAS J214042.81-h581937.4 as YSO 
candidates. 
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3.5. Optical and NIR variability 

Out of the 115 sources (110 MIR + 5 Ha emission sources) mentioned in Sees. 13.21 13.41 and 
we detect a total of 39 sources in B, 46 sources in V, 51 source in R and 63 source in /bands. 
The B VRI magnitudes of all the sources detected in V band are presented in Table [3l We do not 
detect the optical counterparts of the 13 PMS candidates listed in Table [2l Photometric variability, 
a very common characteristic of YSOs, can happen due to various physical processes viz., variation 
in stellar spots or variable obscuration which affect the emission from its photosphere. Since we 
have i? V7 observations separated by several months, we have also looked for the variability in the 
-BV7 bands for all the sources. The median variability in 5 V7 bands are ~ 0.04 mag which is also 
indicative of the accuracy of the optical photometric measurements. We consider a source to be 
variable if the magnitude difference is > 0.1 mag (Table [3]) at the two observing epochs. MIR-5, 
MIR-32, MIR-43, MIR-45 and MIR-49 show variation in 5 V7 bands. In particular MIR-5 with 
WA(Ha) ~ 80 A shows ~ 1 mag variability in B over a timescale of month. 

YSOs also show variation in NIR bands which may or may not have any correlation with 
the optical variation (JEiroa et al.ll2002l ). Both the stellar photosphere and the circumstellar disk 
contribute to the NIR fluxes of PMS stars. Hence, NIR variability can be either due to the 
same physical processes responsible for optical variability or any structural variability of disk and 
circumstellar materials. We compare the NIR mag nitudes from 2MASS (jCutri et al.l |2003| ) and 
from the more recently obtained JHK magnitudes by iBeltran et al.l (|2009l ) of the sources reported 
in both the datasets to identify NIR variability. Sources which have shown variation of ~ 0.1 mag 
in any of the NIR band are listed in the last column of Table [Ij 



3.6. Medium resolution spectra of YSOs 



The Ha emission though a very prominent feature in the spectrum of the YSO is not sufficient 
to identify YSOs because there also exist other kinds of Ha emission line stars viz., Be, Mira stars 
etc. Medium resolution spectroscopy is helpful to distinguish the YSOs from the Be and Mira stars. 
Figure H] presents sample spectra observed with HFOSC of a few YSOs. The prominent features 
e.g. Ha, Call NIR triplet etc. are marked in the respective spectra. We have detected strong 
Ha emission in most of the Class II objects except BRC 38 16 and therefore we do not consider 
it as a YSO in the remainder of the paper. Log of spectroscopic observations and Ha equivalent 
widths (Wa) of the individual YSOs are given in Table [1 BRC 38 6 (MIR-31) is the brightest 
Class II object detected in optical band. From the multiepoch spectroscopic observations we find 
that MIR-31 shows variable Ha emission (Table U}. We have also detected 01 A7773 ( blend of 
01 A7771.96, A7774.18, A7775.40 ), 01 A8446 (blended with Pal8 A8437) Call infrared triplet i.e 
A8498 (blended with Pal6 A8502), A8542 (blended with Pal5 A8545) and A8662 (blended with Pal3 
A8665) and several other Paschen lines e.g. PalO A 9014, Pal l A8862, Pal2 A8750, P al4 A8598, 
Pal7 A8467, Pal9 A8413, Pa20 A8392 in MIR-31. According to lMuzerolle et al.1 (|l998bl ). 01 A7773 
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is a very sensitive indicator of infall in the context of magnetospheric accretion scenario. Lil A6708 
(Wa ~0.7 A) is also present in the spectra of MIR-31. All these spectroscopic features indicate that 
MIR-31 is indeed a YSO which is going through an active accretion phase. We have detected [01] 
A6300 and [01] A6363 and [SII] A6717 and [SII] A6731 emission lines in MIR-45 which are well known 
signatures of outflows from YSOs. We also find a positive correlation between Ha emission and 
MIR excess for YSOs with WA(Ha) < 40 A. We have detected sporadic Ha emission in transitional 
disk object MIR-29 on 3rd November 2008. We observed the object in 2009 also but we do not 
find any further Ha emission. We have detected weak Ha emission from two more transitional disk 
objects e.g. MIR-43 ([GFG2007] 62) and MIR- 76 ([GFG2007] 81). We also obtained the spectra 
of few more transitional disk objects e.g. MIR-1, MIR-2, MIR-9, MIR-18 and MIR-108 and do not 
find Ha emission in any of these objects. 



3.7. Spectral Classification 



We have classified the observed optical spectra by manual inspection of the spectra of target 
stars and also by comparing it wit h the spectra of known spectral ty p e. W e obs erved several stars 



with s pectral type of F to M from Ijacoby et al.l (| 19841 ) , IValdes et al.l (120041 ) and iKirkpatrick et al 



(|l99ll ) with HFOSC using the same set up as the target stars to minimize the effect of instrumental 
signatures and resolution. From preliminary analysis we find that the YSOs are of spectral type K 
and M. 

For more accurate classi fication in the range of spectral type K-M, we use the calibration file 
provided in the SPTCLASS (JHernandez et al.ll2004l ) code for various absorption features viz., Tio 
bands etc. The detailed information of the classification scheme can be found in the SPTCLASS 
codq£|. This scheme measures the equivalent widths or "indices" of the absorption features which 
are sensitive to Tg// using the adjacent blue and red continuum bands and calibrate the indices 
against the spectral type of the stars. Spectral type of the target star is determined from the 
weighted mean of the spectral types obtained from different absorption features of the target star 
excluding the apparently anomalous values. As this classification scheme uses the local continuum 
of the absorption feature, it is insensitive to reddening and the S/N ratio. Further, the spectral 
features which are used for classification are also insensitive to the luminosity class. The errors in 
the spectral type determination come from the measurement of indices and the standard deviation 
of the spectral types obtained from different absorption features. To check the applicability of the 
adopted calibration on HFOSC spectra we classify the HFOSC spectra of a few K and M type 
stars and retrieve the spectral type within ~0.5 subtype. We then use the same calibration scheme 
to classify the YSO sp ectra. We further exp lore the Fe I A6495 A line to classify the spectra 
of early K type stars (JGrav fc Corballvl l2009l ) . The "index" or equivalent width of A6495 line 
was calculated for spectral types F-K, using the HFOSC spectra of stars with known spectral 



'^http : //www . astro . Isa . umich . edu/~ hernandj /SPTclass/sptclass . html | 
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type taken from iJacobv et al.1 (| 19841 ): IValdes et alj (|2004l ). The results of our classification are 
summarized in Table HI W e also tabulate the effective temperatures of YSOs using the table of 
Kenvon &: HartmannI (jl995l ). Using the V-I color of a particular spectral type and our optical 
B VRI observations we ca lculated the A^ assuming the normal interstellar extinction law as given 
bv lBessell &: BrettI (|l988l ). The average extinction towards the globule is estimated to be A^= 2.5 



mag. 



3.8. Rate of accretion from Ha emission 



The Ha emission along with the other Balmer lines, have always been the iconic spectral diag- 
nostics of accretion of circumstellar materials in the YSOs. Magnetospheric accretion models have 
provided satisfactory explanation of the observed spectral cha racteristics of accretion Class ical T 
Tauri Stars (CTTS) e.g. broad asymmetric line profiles etc. (JMuzerolle et al.lll998al . l200ll ). Ac- 



cording to the magnetospheric accretion model the stellar magnetic field disrupts the circumstellar 
disk a t several stella r radii and accreting material fall to the stellar surface along the magnetic field 



lines (|Koenigllll99ll ). Magnetospheric accretion models suggest that the hydrogen emission lines 
form in the infall zone, so that the blue shifted asymmetric emission line profiles arise due to partial 
obscuration of the flow by the inner part of the accretion disc and the red shifted absorption pro- 



files re sult from infalling material at near free-fall velocities on the stellar surface (JMuzerolle et al 



1998al ). Strength of the Ha emission line has been widely used to distinguish the actively accret- 
ing CTTS frorn their less active and evolved counterparts i.e Weak-lined T Tauri Stars (WTTS). 
White fc Basril j^Ooi) proposed to classify a YSO as CTTS if Wa (Ha) > 3 A for K0-K5 stars, 
Wa (Ha) > 10 A for K7-M2.5 stars, Wa (Ha) > 20 A for M3-M5.5 stars, and Wa (Ha) > 40 A for 
M6-M7.5 stars. They also proposed that the stars with full width of Ha emission line at 10 % of 
the peak intensity (Ha[10%]) greater than 270 km s~^ are accreting CTTS independent of spectral 
type. We i dentify 9 sta r s as C TTS and 5 stars as WTTS according to the above mentioned criteria 
(Table El). iNatta et al.1 (I2OO4I ) showed that Ha[10%] > 200 km s~i can be used to distinguish the 
accretors and non accretors as well as to derive an approximate estimate of the accretion rate M^c 
quantitatively by the following relation, logM^c = -12.89 (±0.3) + 9.7 (±0.7) x 10"^ Ha[10%] 
where Ha[10%] is in km s~^ and M^c is in M0yr~^. In order to estimate the accretion rates for some 
of the well resolved Ha line profile we fitted Gaussian profiles to the Ha emission lines and calculate 
the equivalent width and FWHM of each line profile. The measured FWHM of each profile has 
been deconvolved assuming a Gaussian instrumental profile. We calculate the Ha [10%] for these 
objects and transformed them into velocity units. We did not correct for the underlying photo- 
spheric absorption in Ha if there is any. Table |4] presents the accretion rates M^c, calculated using 
the above formula along with the measured values of Ha[10%]. Based on the observed multiepoch 
spectra we find variable accretion rates in BRC 38 1, MIR-5, MIR-31, MIR-32 and MIR-45. 
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3.9. Optical Color- Magnitude Diagram 



FigureOpresents the extinction corrected V vs. V-I diagram of tiie YSOs near SFO 38. We have 
corrected for the extinction of individual objects with known spectral types and used an average 
value of Av=2.5 mag when the spectral t ype is not known. In order to estimate the mass and age of 
the Y SOs we overplot the isochrones from lSiess et al.l (2000) assuming a distance of 750 pc lMatthews 
(|l979l ). We find that most of the YSOs (Class II with(out) Ha) have ages between 1-5 Myr and 
masses in the range 0.3-2.2 M.q. We note that several transitional disk objects are located to the 
left of the ZAMS line. A closer examination of the position of the transitional disk objects in the 
MIR color-color diagram (Fig. [3j) reveals that 6 out of 7 transitional disk objects which are situated 
to the left of ZAMS line have [5.8]-[8] color excess < 0.5 and 4 out of 5 transitional disk objects 
which are situated to the right of ZAMS line have [5.8]-[8] color excess > 0.5. The transitional disk 
objects to the left of the ZAMS line are more like Class III objects with smaller values of color 
excess. We have used an average A^=2.5 mag for extinction correction for these objects because 
the spectral types are not known. There are several factors that could lead the transitional disk 
objects to the anomalous position in Fig. [5] e.g. extinction correction, uncertainties in distance and 
PMS evolution model used. An extinction correction with a smaller A^ (<1 mag) brings all the 
transitional disk objects closer to, or t o the right of the ZAMS line. The distance of IC 1396 is in the 
range of 615 pc (jde Zeeuw et al.lll999l ) to 900 pc (jContreras et al.ll2002l ) and we adopt a distance of 
750 pc for SFO 38 to make our results comparable to most published results. This uncertainty in 
the distance corresponds to an uncertainty of it 0.4 mag in the absolute magnitude. Adoption of 
a distance of 900 pc brings up the transitional disk objects somewhat closer to the ZAMS line but 
these objects still occupy positions below the ZAMS line. While different theoretical models differ 
in the details of the PMS evolution, close to the ZAMS line all these models agree with each other 
and with the empirical ZAMS of the young star clusters. Therefore, the anomalous positions of the 
transitional disk objects with respect to the ZAMS track in Fig.[5]is primarily due to overcorrection 
of the extinction. Estimation of m ass and age from the colo r-magnitude diagram (CMD) is strongly 
dependent on the adopted model (JHillenbrand et al.l 120081 ). The age s can differ by as much as 10 
Myr or more when two different models are used on the same dataset (llkeda et al.ll2008l ). The main 
source of error in the V vs. V-I diagram is the estimation of extinction towards the globule. In the 
last column of Table H] we list the ratios of E(V-I)/E(B-V) which are different from the standard 
value of 1.6 and indicate that normal extinction law may not hold good for star forming regions 
due to plenty of patchy obscuring materials. Although the NIR JHK bands suffer comparatively 
less extinction than the optical BVRI hand but for the YSOs, thermal emission also contributes 
to JHK bands resulting in excess fluxes compared to the main sequence counterparts without disk. 
For a given model we consider the derived values of the masses and ages as representative values 
and the age spread of YSOs as more reliable estimate. 
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3.10. SED of YSOs with MIPS 24 /xm detection 

In order to interpret the observed SEDs for the YSOs detected at 24 fim and better charac- 
terize these sources we have explored the archive of two-dimensional (2D) axisymmetric radiative 
transfer models of protostars calculated for a large range of prot ostellar masses, accretion rates, disk 



masses and disk orientations created by iRobitaille et alJ (|2007l ). This archive also provides a linear 



regression tool which can select all model SEDs that fit the observed SED better than a specified 
X^. Each SED is characterized by a set of model parameters, such as stellar mass, temperature, 
and age, envelope accretion rate, disk mass, and envelope inner radius. We have used this online 
tool to generate models, which fit the observed SEDs for the 14 YSOs which have been detected 
at 24 /xm. We restricted the SED fitting tool to explore only distances between 650 and 850 pc. 
Below we present and briefly discuss the best fit models. 

Figure [6] shows results of detailed modeling of the observed SED in the mid-infrared (and 
sub-mm for MIR-50) for all the YSOs which were detected at 24 ^um. A major criticism against 
the use of these models has been the non-uniqueness of the solutions obtained from the model 
library. However we find that for all sources the best fit models are distinctively better than the 
next few models in reproducing all the observed fiux densities. Table [S] presents the parameters 
corresponding to the best fit models for the fourteen sources which we fitted with the accretion 
disk models. 



Robitaille et al.l (120061 ) presented a classification scheme which is essentially analogous to the 
Class scheme, but refers to the actual evolutionary stage of the object based on its physical prop- 
erties like disk r nass and envelop e accretion rates rather than the slope of its near /mid- IR SED. 



According to the IRobitaille et al.l classification scheme, MIR-50 and MIR-54, are stage O/I objects, 
while all the others are stage II objects. However, the YSOs MIR-34, 36, 48, 55 and 59 were iden- 
tifled as Class O/I objects based on the color-color diagrams. Next we consider the fact that except 
for MIR-34 all the sources have rather high visual extinction (^v >40) and the MIR color-color 
diagrams we had used to c lassify the YSOs did n ot correct for the reddening. Using standard 



interstellar extinction laws (jlndebetouw et al.ll2005l )and the A^ derived from the models, we flnd 



that when de-reddened, the sources MIR-48, 55 and 59 occupy the region typically occupied by 
Class II objects. However for MIR-34 and 36 the difference in the inferred evolutionary stage can 
not be explained only in terms of reddening. 

Based on the results of fltting the observed SEDs with axially symmetric accretion models we 
find that all the YSOs analyzed here have low to intermediate-mass central stars. Figure [2] shows 
that consistent with our expectations the sources located within the globule have ^v>20 mag, while 
the sources lying at the periphery of the bright rim have A^ between 3-5 mag. In Sec. 13.111 we 
discuss the two most luminous Class O/I sources MIR-50 and MIR-54 in detail. 

The ages estimated for these fourteen YSOs from the SED models and the location of these 
sources in BRC 38 show that most of the sources closer to the bright rim have ages between 5- 
8 Myr. The sources within the BRC and lying almost along the axis of the BRC have ages between 
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0.1-1 Myr, with the protostehar cluster at the position of IRAS 21391+5802 being the youngest 
with an age of -^0.1 Myr. The age spread estimated for the YSOs (mostly Class II and above) from 
optical CMD diagram is of the range 1-8 Myr. However, the ages estimated from both methods 
are higher than the ages obtained from the NIR CMD by iGetman et al.l (J2007l ) . 



3.11. The protostellar cluster at IRAS 21391+5802 

Interferometric millimeter continuum images have shown that the most luminous source IRAS 
21391+5802 in the IC 1396 region can be resolved into three components, named as BIMA 1, 
BIMA 2 and BIMA 3 (JBeltran et al.ll2002l ). Radio continuum images at 3.6 cm observed with 
VLA have also iden tified three peaks VLA 1-3, coincident with the millimeter continuum peaks 
(JBeltran et al.l |2002| ) . Of these three sources VLA 3 is the strongest emitter of radio continuum. 
Recent sub-arc second reso l ution millimeter observations have identified three components within 
BIMA 2 itself iNeri et al.l ([2003). The SCUBA images at 450 and 8 50 Mm do not reso lve the 
components of continuum emission associated with IRAS 21391+5802 (JMorgan et al.l 120081 ). 



We identify MIR-50 and 54 as the mid-infrared counterparts of BIMA 2 and BIMA 3 and 
do not detect any source associated with BIMA 1. The source MIR-36 is located 3'.'8 away from 
BIMA 1 and hence is most likely not associated with BIMA 1. We have used the continuum 
fluxes available in literature to construct the observed SEDs for MIR-50 and MIR-54. In order 
to make a realistic estimate of the flux densities at 450 and 850 ^m from each of BIMA 2 and 
BIMA 3, we have split up the observed total flux densities at these wavelengths a ccording to the 
ratio of flux densities of these sources measured at 1.2 mm by lBeltran et al. (120021). Fitting of the 
observed flux densities with the axially symmetric accretion-based models (jRobitaille et al.l 120071 ) 
shows that MIR-50 (BIMA 2) is a protostellar object with L = 197 Lq and a mass of 5.97 M©. Both 
mass and luminosity values match reasonably well with the literature. We note however that the 
model fitted does not reproduce the interferometric fl ux densities at 1.2 and 3.1 mm well. This is 
consistent with the conclusions of lBeltran et al.l (|2002l ) that the interferometric observations miss a 
significant amount of emi ssion. The best fit mo del for MIR-54 corresponds to a mass of 1.5 Mq and 
a luminosity of 33.4 Lq. ICodella et al.l (J200l|) had estimated an integrated luminosity of 233 Lq 
including both BIMA 2 and BIMA 3 and here we derive a total luminosity of 230 Lq. Based on 
the SED fitting performed in this paper, MIR-54 or BIMA 3 is a lo w mass Cla s s O/I protostar and 
MIR-50 or BIMA 2 is a Class O/I intermediate-mass star, although iNeri et al.l ((20071) have already 
shown BIMA 2 as having three components. The SED models derive the age of both these sources 
to be ~ 10^ yr. 
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Discussion: RDI mechanism at work in SFO 38 



Based on the observational evidences in the form of YSOs, an outflow-driving IRAS source 
and several HH objects, SFO 38 is unambiguously an active star forming region. There is however 
considerable disagreement regarding whether the star formation is triggered or influenced by the 
nearby OB association. The two m ost recent studies of S FO 38 based on Chandra X-r ay and Syitzer 



Mid- IR observations (up to 5.8 jum lGetman et al.ll2007l ) and deep NIR observations (JBeltran et al 



20091 ) arrive at diametrically opposite conclusions. Considering the spatial and temporal distribu- 
ti on of the fewer YSO s selected by using MID-IR color-color diagram ( [3.6]-[4.5] vs. [4.5]-[5.8] 
) , iGetman et al.l (|2007l) argued in favor of t riggered star formation by Radiation Driven Implosion 
(RDI) model. However iBeltran et al.l (|2009l ) did not find either (a) enhanced star formation activity 
towards the southern rim of SFO 38 or (h) systematic gradient in the evolutionary stage_oftheYSOs 
using the deep NIR imaging. Based on the absence of NIR excess from the YSOs iBeltran et al 



(|2009l ) concluded that the sources closer to the O type star appear to be more evolved primarily due 



to the modification of their environme nt by the inten se UV field. Thus, according to IBeltran et al 



the evolutionary sequence observe d bvlGetman et al.l should not necessarily be interpreted as due to 
triggered star formation. However IBeltran et al.l also concede that NIR-only studies are inadequate 
to identify all the YSOs in the region. 

We explore two possible scenarios which could have been the dominant driving mechanism to 
initiate the formation of stars in SFO 38 viz., coeval star formation throughout the H ii region and 
triggered star formation by the RDI mechanism. The assumption that the stars near the globule 
were formed simultaneously with the other stars of the OB association suggests that the distribution 
of YSOs should be correlated with the distribution pattern of the other massive to low mass stars 
of the OB association. The observed elongated distribution of YSOs along two different axes in 
SFO 38 is quite different from t he distribution of the known members of Cep 0B2 association 
(Fig. 22 of lde Zeeuw et al.lll999l ). Moreover the ages of the outflow driving IRAS sou rce and the 



associated protostellar cluster are significantly less than the dynamical age of 3 Myr (jPatel et al. 



19951 ) of the H ii region created by the massive OB type stars. We thus rule out the possibility 
of coeval star formation in SFO 38 with the members of Cep 0B2 association, although some of 
the Ha emitting sources could have been formed prior to the arrival of the shock front. Next, we 
investigate the other possibility, based on the RDI model, by comparing our observations with the 
predicted signatures of the model. The RDI model suggests enhancement of density along the axis 
of the UV irradiated globule and an age sequence of YSOs with the oldest YSO lying closer to the 
OB association. 

Sequential st ar formation i s cons idered to be one of the key signatures of triggered star for- 
mation via RDI. IBeltran et al.l (J2009l ) have cautioned against the use of the spatial distribution 
of YSOs as a proof of triggered star formation, by arguing that intense UV fields close enough 
to the OB associations tend to destroy the protoplanetary disks around the YSOs and thus make 
them appear to be older. However we find this argument unt enable for the case of SFO 38, since 



the typical radius of infiuence of O stars is found to be < 1 pc (JBalog et al.l 120071 : iHernandez et al 
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20081 : iMercer et al.ll2009l ). while the mean projected distance of the nearest O and B -type stars 
close to SFO 38 is ~9 pc. Thus the spatial distribution of the YSOs discussed in SFO 38 signify 
the temporal evolution of the YSOs and not due to the destruction of protoplanetary atmospheres 
under the influence of the UV radiation from the OB association. 

The bright-rimmed cloud SFO 37, was found t o have an an exem plary bow shaped geometry 



and show a clean evolutionary sequence of YSOs (jlkeda et al.ll2008l ). SFO 38 in contrast has a 
more complicated geometry. As noted earlier the Ha emission from SFO 38 shows an asymmetry 
relative to the axis connecting the globule to the nearest O star HD 206267. Further, in SFO 38 
though the YSOs clearly show an age segregation, they are not distributed along a single straight 
line. It is however possible to identify another direction along which some of the YSOs are most 
likely aligned. The asymmetry of Ha emission and the non-linearity of the distribution of YSOs 
suggest that the evolution of SFO 38 is influenced by the nearby B star HD 206773 as well. From 
the Ha line image (Figure [2|) HD 206773 seems to be ionizing the eastern part of the rim more 
efficiently than the extent to wh ich HD 206267 ha s ionized the western part of the rim. Using the 



13/ 



CO J = 1-0 map of SFO 38, iPatel et al.l (119951 ) also found that the eastern and western wings 



of SFO 38 differ in velocity by about 1.5 km s ^ relative to the velocity of the dense core of the 
globule. This suggests that HD 206773 plays an active role in the present dynamics of the globule. 

We have identified a spatio-temporal gradient along the directions of both HD 206267 and 
HD 206773. However the spatial gradient is not continuous, there is a distinguishable gap between 
the two classes of YSOs e.g. the projected distance of Class I/II sources from the approximate 
center of the distribution of Class I sources is ~ 0.3 pc and this gap can be interpreted as due to 
the difference in the arrival times of the shock within the framework of the RDI model. If we consider 
the age spread to be ~ a few Myr (comparing the estimated ages from SED fitting and optical 
CMD) and the spatial separation between the MIR-50 (Class O/I) and MIR-29 (Class II) to be ~ 
0.6 pc, we obtain the speed of shock propagation to be 0.1-0.3 km s~^. The esti mated shock spee d 



is consistent with the shock velocities obtained from 3-D numerical simulations (|Miao et al.ll2006l ). 
The estimated speed of shock propagation are approximate since we use only projected distances 
between sources and do not consider the time lapse between the arrival of the shockw ave and the 



trigg ering of star formation in the globule, which from simulations is ~ 0.375 Myr (JMiao et al 



20061 ^ 



Thus, our observations and subsequent estimate of the evolutionary stages of the YSOs embed- 
ded in SFO 38, are consistent with a scenario of triggered star formation via the RDI mechanism 
powered by two different ionizing sources, HD 206267 and HD 206773. 



5. Summary 

We have presented multiwavelength photometric study of SFO 38 from optical BVRI to Spitzer 
IRAC and MIPS observations along with optical spectroscopy of the selected objects. We summarize 
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our work as follows. 

• A total of 40 YSOs (Classes O/I/II) and 13 YSO candidates are identified based on MIR color 
indices and we also confirm Ha emission from 2 YSOs which are not known in the literature. 
We further identify 4 additional YSOs based on Ha emission and thus we find 44 YSOs in 
and around SFO 38. 

• Ha emission line YSOs are mostly K-M type stars with an age spread of 1-8 Myr and mass 
range of .3-2.2 M0Some of the YSOs show photometric variation in optical and NIR bands and 
also variable Ha emission in the medium resolution spectra. Mass accretion rates estimated 
form broad Ha line profiles are of the order of 10^*^ to 10 ^^^ Mq yr^^. 

• Mass, luminosity and age for the different components of the protostellar cluster at IRAS 
21391+5802 were derived by SED fitting. 

• Continuum subtracted Ha line image show asymmetric Ha emission at the bright-rim. Two 
OB type stars e.g. HD 206267 (06.5 ) and HD 206773 (BOV) are proposed as the potential 
ionizing sources. 

• Class II to Class O/I objects are distributed from the rim to the core part of the globule 
respectively. It is possible to identify at least two different axes of elongation of the YSO 
distribution. The spatial gap between the two different classes of YSOs is consistent with 
the difference in the arrival time of a shockwave propapagating into the globule. The spatio- 
temporal gradient in the distribution of YSOs along two different axes that are parallel to 
one of the either ionizing star indicate triggered star formation due to the Radiation Driven 
Implosion model. 

This work is based [in part] on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract 
with NASA. This research has made use of the SIMBAD data base, operated at CDS, Strasbourg, 
France. This publication makes use of the data products from the 2MASS, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronautics and Space Administration and the 
National Science Foundation. RC would like to thank Dr. P. S. Parihar and Ramya S. for useful 
discussions. 
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Fig. 1. — IRAC-MIPS color-composite image of SFO 38 using 3.6 ^m (Blue), 8 fim (Green) and 



24 ^m (Red). The image is centered at a2ooo 
5'x5' (a x6). 
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Fig. 2. — Continuum subtracted Ha emission line image of SFO 38. Class O/I to Class II YSOs are 
overplotted with different symbols. The directions towards HD 206267 and HD 206773 are marked 
with two different arrows. 
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Fig. 3. — Color-colo r diag r ams for all MIR so urces . Approximate c lassifi cation criteria adopted 
from iMegeath et al.l (J2004l ) , iFang et al.l (120091 ) and iMuzerolle et al.l (120041 ) are shown in the two 
panels, (a) MIR color-color diagram based on Spitzer 3-band IRAC and MIPS ph otometry and (b) 



4_band IRAC color-color diagram. Dashed lines in the two panels are taken from iHartmann et al 



(120051). Reddeniiig vect or corresponding to the extinction laws given by the fitted function from 
Indebetouw et al.l (|20od ) is shown in (b). The dashed lines [3.6]-[4.5] = 0.7, [4.5]-[5.8] = 0.7 and 
[5.8]-[8.0] = 0.7 discriminate Class II sources from Class I/O sources, and [3.6]-[4.5] = 0.2, [4.5]- 
[5.8] = 0.2 and [5.8]-[8.0] = 0.2 discriminate Class III from Class II and Transitional Disk sources, 
(c) NIR an d Spitzer MI R color-color diagram. Solid line is the reddening vector due to interstellar 
extinction (JTapialll98ll ). Sources to the right of the reddening vector are YSOs with NIR excess. 
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Fig. 4. — Sample spectra of YSOs in SFO 38 obtained with HFOSC instrument. The strong 
features e.g. Ha emission (A6563A ), Call infrared triplet emission , 01 line at A7773 A and [01] 
line at A6300 A are marked in the spectra. Atmospheric features are indicated with ©. Two spectra 
of MIR-29 taken at two different epochs show variable Ha emission (see text for details). 
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Fig. 5. — Extinction corrected V vs. V-I Color-Magnitude diagram of YSOs near S FO 38. Ha 
emissi on YSOs are shown by filled circles. The ZAMS and 0.6, 1, 5, 10 Myr isochrones of lSiess et ah 
(J2000l ) are shown with the evolutionary tracks for masses from 0.3 to 2.2 M©, at an adopted distance 
of 750 pc. 
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Fig. 6. — SED fits for the PMS objects using axisymmetric radiation transfer models. The fihed 
circles indicate the measured fluxes and uncertainties. The filled triangles correspond to the upper 
limits of flux at these wavelengths. The continuous black line represents the best-fitting SED, and 
the dashed line shows the stellar photosphere corresponding to the central source of the best fitting 
model, in the absence of circumstellar dust (but including interstellar extinction). 



Table 1. Results of Near- & Mid-IR photometry of the Mid-infrared sources in SFO 38 



Source 


tl2000 


•52000 


J 


H 


^8 


^3.6 


^4.5 


^5.8 


Fa.a 


^24 




Class 




Commc 


nts 


SF038 












mjy 


mJy 


mJy 


m,Ty 


m 


Jy 










MIR-1 


21:40:24.97 


58:16:10.5 


14.48±0.04 


13.72± 0.05 


13.42± 0.05 


1.67±0.01 


l.lOiO.Ol 


().96±0.()2 


1.03±0.02 






Trans. 


JHKs" 




MIR-2 


21:40:25.39 


58 


15:47.3 


15.35±0.06 


14.51± 0.07 


14.15± 0.07 


0.94±0.01 


0.63±0.01 


0.68±0.02 


0.98±0.02 






Trans. 








MIR-3 


21:40:25.84 


58 


15:26.1 


15.62±0.07 


14.80± 0.07 


14.62± 0.10 


0.52±0.01 


0.31±0.00 


0.23±0.02 


0.07±0.02 






III/S 








MIR-4 


21:40:26.84 


58 


14:51.8 


16.26±0.10 


15.67± 0.15 


15.06± 0.15 


0.32±0.00 


0.20±0.00 


0.17±0.02 


0.11±0.02 






III/S 








MIR-5 


21:40:27.29 


58:14:21.14 


14.30±0.04 


13.30± 0.04 


12.88± 0.04 


6.22±0.()5 












II 


Ha 


X-ray 




MIR-6 


21:40:27.38 


58 


16:21.3 








0.07±0.00 


0.03±0.00 


0.51±0.02 


2.00±0.03 






I/II 








MIR- 7 


21:40:27.97 


58 


15:14.1 


14.51±0.04 


13.41± 0.04 


12.94± 0.04 


3.89±0.02 


3.52±0.02 


3.44±0.03 


4.33±0.03 






II 


Ha 


J 




MIR-8 


21:40:28.35 


58 


16:42.3 


11.59±0.03 


10.78± 0.03 


10.43± 0.02 


23.62±0.07 


14.18±0.04 


10.94±0.06 


6.63±0.04 






III/S 








MIR- 9 


21:40:30.63 


58 


15:00.0 


15.31±0.05 


14.90± 0.07 


14.70± 0.11 


0.49±0.01 


0.32±0.01 


0.28±0.01 


0.22±0.02 






Trans. 








MIR- 10 


21:40:31.00 


58 


15:09.2 


13.44±0.03 


12.97± 0.04 


12.79± 0.04 


2.61±0.02 


1.63±0.01 


1.29±0.02 








III/S 








MIR- 11 


21:40:31.56 


58 


17:55.1 


14.03±0.03 


12.89± 0.03 


12.39± 0.03 


6.93±0.03 


6.93±0.03 


9.15±0.06 


15.17±0.12 


20.46±0 


08 


II 


Ha 


X-ray, 


JHKe 


MIR-12 


21:40:31.97 


58 


14:56.0 


16.72±0.16 


15.78± 0.15 


15.41± 0.21 


0.41±0.00 


0.28±0.00 


0.26±0.01 








III/S 








MIR-13 


21:40:32.23 


58 


16:53.6 


14.13±0.03 


12.44± 0.04 


11.71± 0.03 


10.82±0.04 


6.39±0.02 


6.54±0.05 


7.36±0.05 






Trans. 








MIR- 14 


21:40:32.46 


58 


13:47.2 


13.40±0.03 


12.83± 0.03 


12.63± 0.03 


2.78±0.03 




1.30±0.05 








III/S 








MIR-15 


21:40:32.71 


58 


17:23.5 


12.19±0.03 


11.55± 0.03 


11.36± 0.03 


9.21±0.04 


5.50±0.02 


4.04±0.04 


1.71±0.04 






III/S 








MIR-16'' 


21:40:32.84 


58 


15:47.8 


17.42±0.01 


16.61±0.01 


16.09±0.01 


1.42±0.01 


0.73±0.01 


9.50±0.05 


32.10±0.11 






I/II 








MIR-17 


21:40:33.64 


58 


14:44.6 


16.54±0.14 


15.81± 0.18 


15.07± 0.00 


0.21±0.00 


0.13±0.00 


O.lliO.Ol 


0.16±0.02 






Galaxy 








MIR-18 


21:40:33.69 


58 


14:59.5 


13.06±0.03 


12.09± 0.03 


11.75± 0.02 


8.13±0.03 


4.82±0.02 


4.72±0.04 


5.60±0.03 






Trans. 








MIR-19 


21:40:34.05 


58 


18:08.1 


15.70±0.06 


13.40± 0.04 


12.30± 0.02 


7.41±().()3 


5.16±0.02 


4.95±0.04 


4.03±0.04 






Trans. 






1 


MIR-20 


21:40:34.30 


58 


18:23.6 








0.31±().()() 


().S4±().()1 


1.03±0.02 


1.24±0.04 






H2 






KS 


MIR-21 


21:40:34.59 


58 


16:11.4 








0.08±().()0 


0.05±().0() 


1.33±0.03 








III/S 






C5 


MIR-22'' 


21:40:34.76 


58 


15:20.9 


17.70±0.01 


16.68±0.01 


16.04±0.01 


0.25±().()1 


().34±().()1 










III/S 






1 


MIR-23 


21:40:34.79 


58 


16:10.8 














4.47±0.05 














MIR-24 


21:40:35.02 


58 


18:22.2 


16.48±. . . 


15.72±. . . 


14.90± 0.21 


0.62±0.01 


2.68±0.01 


1.84±0.03 


1.42±0.03 






H2 








MIR-25'' 


21:40:35.11 


58 


15:30.9 


18.21±0.02 


16.S5±0.01 


16.16±0.02 


0.92±0.01 


0.48±0.01 


5.32±0.04 


18.26±0.08 






I/II 








MIR-26 


21:40:35.20 


58 


14:05.0 


13.11±0.03 


12.61± 0.03 


12.45± 0.03 


3.55±0.02 


2.40±0.01 


1.59±0.02 


0.84±0.02 






III/S 








MIR-27 


21:40:35.68 


58 


18:21.3 


16.55±... 


15.73±... 


14.79± 0.17 


0.77±0.01 


3.61±0.02 


2.20±0.03 


1.91±0.03 






H2 








MIR-28 


21:40:36.46 


58 


15:23.6 


15.97±0.10 


14.73± 0.08 


14.04± 0.07 


2.27±0.01 


1.43±0.01 


6.68±0.04 


19.47±0.08 






I/II 








MIR-29 


21:40:36.52 


58 


13:45.8 


13.51±0.02 


12.58± 0.03 


12.24± 0.03 


4.69±0.02 


3.14±0.01 


2.17±0.03 


2.64±0.03 


44.29±0 


09 


Trans. 


Ha 


X-ray 




MIR-30 


21:40:36.65 


58 


18:23.0 


17.12±. . . 


15.18±0.10 


13.93±0.07 


1.84±0.01 


1.55±0.01 


1.42±0.03 


1.21±0.03 






II 








MIR-31 


21:40:36.88 


58 


14:37.8 


11.90±0.02 


10.89± 0.03 


10.23± 0.02 


54.30±0.14 


48.79±0.10 


38.25±0.12 


34.45±0.10 


39.36±0 


07 


II 


Ha 


X-ray 




MIR-32 


21:40:37.02 


58 


15:02.9 


14.27±0.03 


13.28± 0.04 


12.82± 0.03 


4.01±0.02 


3.07±0.01 


2.86±0.03 


3.87±0.04 






II 


Ha 


Ks 




MIR-33 


21:40:37.16 


58 


15:46.6 


17.01±. . . 


15.80±0.16 


14.55±0.10 


1.33±0.01 


1.28±0.01 


0.71±0.03 








III/S 








MIR-34'' 


21:40:37.54 


58 


15:39.3 


18.07±0.02 


16.98±0.01 


16.24±0.02 


0.49±0.01 


0.77±0.01 


1.05±0.03 


1.63±0.04 


17.15±0 


06 


O/I 








MIR-35 


21:40:39.05 


58 


16:03.3 








0.41±0.01 


1.04±0.01 


2.19±0.03 








H2 








MIR-36'* 


21:40:39.60 


58 


16:09.0 




>17.6 


15.14±0.16 


10.47±0.03 


36.11±0.08 


59.09±0.17 


82.89±0.22 


239.10±0 


10 


O/I 


X-r 


ay 




MIR-37 


21:40:39.80 


58 


18:24.9 








0.53±0.01 


0.76±0.01 


1.15±0.02 


1.39±0.03 






O/I 








MIR-38 


21:40:39.84 


58 


18:34.8 


18.29±. . . 


15.30±0.12 


13.36±0.04 


6.88±0.03 


8.19±0.03 


9.95±0.06 


13.49±0.06 


30.69±0 


07 


II 


X-r 


ay 




MIR-39 


21:40:40.15 


58 


16:18.2 








1.36±0.01 


2.65±0.01 


3.24±0.04 


4.25±0.05 






O/I 








MIR-40 


21:40:40.43 


58 


18:26.1 


18.28±. . . 


15.64±0.15 


13.85 ±0.06 


2.35±0.01 


2.05±0.01 


1.64±0.03 


0.76±0.02 






III/S 








MIR-41 


21:40:40.96 


58 


16:12.0 


17.55±... 


15.89±. . . 


13.59± 0.06 


13.04±0.04 


23.41±0.06 


29.43±0.10 


42.25±0.15 






O/I 








MIR-42 


21:40:40.98 


58 


16:52.4 


17.93±... 


15.34±0.14 


13.93±0.06 


1.41±0.01 


1.04±0.01 


0.83±0.02 


0.37±0.03 






III/S 








MIR-43 


21:40:41.07 


58 


13:58.9 


12.96±0.02 


12.08± 0.03 


11.77± 0.02 


6.80±0.02 


4.58±0.02 


3.19±0.03 


1.84±0.03 






III/S 


Ha 


X-ray 




MIR-44 


21:40:41.10 


58 


17:54.1 


14.20±0.04 


13.46± 0.04 


13.21± 0.03 


1.89±0.01 


1.14±0.01 


0.73±0.02 


0.37±0.02 






III/S 








MIR-45 


21:40:41.14 


58 


15:11.2 


12.97±0.03 


11.61± 0.04 


10.68± 0.02 


66.38±0.17 


65.98±0.14 


81.84±0.21 


124.60±0.31 


194.40±0 


09 


II 


Ha 


X-ray 




MIR-46 


21:40:41.27 


58 


17:14.2 








0.06±0.00 


0.03±0.00 


0.65±0.02 


2.61±0.03 






I/II 








MIR-47 


21:40:41.40 


58 


14:17.5 


15.48±. . . 


15.50±0.12 


14.63±0.00 


0.34±0.00 


0.21±0.00 


0.16±0.01 


0.10±0.02 






III/S 








MIR-48 


21:40:41.43 


58 


16:37.8 


18.30±... 


15.60±0.15 


12.79±0.03 


51.81±0.14 


109.40±0.20 


158.50±0.34 


199.60±0.44 


1362. 00±0 


26 


O/I 








MIR-49 


21:40:41.54 


58 


14:25.5 


13.65±0.03 


12.62± 0.03 


12.17± 0.03 


8.70±0.03 


7.38±0.03 


6.65±0.04 


6.42±0.04 






II 


Ha 


X-ray 





Table 1 — Continued 



Source 


c^2000 


•52000 


J 


H 


Ks 


F3.6 


Fi.5 


^5.8 


fs.o F2i 


Class 


Comments 




SF038 












mjy 


mJy 


mJy 


m,Ty m 


Jy 






MIR-50 


21:40:41.71 


58:16:12.6 








10.84±0.04 


78.01±0.16 


128.00±0.29 


137.00±0.34 1871. 00±0 


31 O/I 


X-ray 




MIR-51 


21:40:41.89 


58 


15:23.0 


15.68±0.09 


14.30± 0.06 


13.65± 0.05 


2.88±0.02 


2.30±0.01 


3.75±0.04 


8.33±0.06 


I/II 


X-ray, JHKs 




MIR-52 


21:40:41.96 


58 


13:49.1 


15.50±0.06 


15.00± 0.09 


14.98± 0.14 


0.40±0.00 


0.27±0.00 


0.19±0.01 


0.10±0.02 


III/S 






MIR-53 


21:40:42.35 


58 


18:40.1 








0.48±0.00 


2.16±0.01 


1.99±0.03 


1.26±0.02 


H2 






MIR-54 


21:40:42.77 


58 


16:01.1 








20.16±0.06 


75.86±0.15 


113.50±0.27 


117.60±0.31 2484.00±0 


38 I 


X-ray 




MIR-SS'' 


21:40:43.32 


58 


17:16.5 




>17.6 


15.31±0.11 


1.65±0.01 


2.23±0.01 


2.72±0.03 


4.31±0.03 145.40±0 


08 O/I 






MIR-56 


21:40:43.34 


58 


18:42.6 


16.32±0.13 


14.42± 0.05 


13.55± 0.04 


1.93±0.01 


1.39±0.01 


0.91±0.02 


0.26±0.03 


III/S 






MIR-57 


21:40:43.62 


58 


16:18.8 


17.89±. . . 


16.09±. . . 


13.51± 0.07 


33.21±0.09 


68.55±0.14 


82.72±0.21 


91.38±0.24 


O/I 


X-ray 




MIR-58 


21:40:43.65 


58 


14:18.4 


15.55±0.08 


15.05± 0.10 


14.89± 0.14 


0.41±0.00 


0.28±0.00 


0.19±0.01 


0.14±0.02 


Trans. 






MIR-59 


21:40:44.23 


58 


16:52.2 








1.36±0.01 


4.43±0.02 


9.13±0.05 


10.52±0.06 100.90±0 


07 O/I 






MIR-60 


21:40:44.26 


58 


14:12.4 


15.57±0.06 


15.00± 0.09 


14.69± 0.12 


0.45±0.00 


0.30±0.00 


0.22±0.01 


0.09±0.02 


III/S 






MIR-61 


21:40:44.31 


58 


15:13.2 


16.05±0.10 


14.59± 0.07 


13.60± 0.05 


3.78±0.02 


4.10±0.02 


7.89±0.05 


15.74±0.08 


I/II 


X-ray, JHKs 




MIR-62 


21:40:44.64 


58 


17:04.9 


17.83±. . . 


15.75±0.15 


14.62±0.10 


2.16±0.01 


2.38±0.01 


2.51±0.03 


3.06±0.03 


II 






MIR-63 


21:40:44.82 


58 


16:05.0 


15.86±. . . 


14.28±0.06 


12.89±0.05 


7.08±0.03 


7.17±0.03 


6.83±0.05 


7.25±0.05 


II 


X-ray, Kg 




MIR-64 


21:40:44.83 


58 


15:03.3 


14.62±0.04 


13.35± 0.04 


12.66± 0.03 


5.64±0.02 


5.23±0.02 


7.28±0.05 


12.94±0.06 


I/II 


Ha, X-ray 




MIR-65 


21:40:44.89 


58 


17:47.7 








0.42±0.00 


2.56±0.01 


2.01±0.02 


2.29±0.03 


H2 






MIR-ee*" 


21:40:45.11 


58 


15:24.0 


18.46±0.02 


16.74±0.01 


15.71±0.01 


0.59±0.01 


0.44±0.01 


1.73±0.03 


5.26±0.05 


I/II 






MIR-67 


21:40:45.16 


58 


15:59.6 








5.73±0.02 


9.26±0.03 


11.46±0.06 


11.78±0.07 


I 


X-ray 




MIR-68'' 


21:40:45.24 


58 


16:41.3 


17.38±0.01 


16.29± 0.01 


15.53±0.01 


0.43±0.01 


0.42±0.01 


0.41±0.02 


0.68±0.02 


Galaxy 




1 


MIR-69 


21:40:45.29 


58 


16:53.3 








0.53±0.01 


0.61±0.01 


0.60±0.02 




III/S 




K3 


MIR-70'' 


21:40:45.41 


58 


16:45.5 


20.38±0.15 


17.48±0.02 


15.96±0.01 


0.33±0.00 


0.27±0.01 


0.48±0.02 


1.14±0.02 


Galaxy 




MIR-71 


21:40:45.48 


58 


15:11.5 


14.65±0.05 


13.71± 0.06 


13.11± 0.04 


3.17±0.02 


2.61±0.01 


3.41±0.04 


5.73±0.05 


I/II 


X-ray, Ks 




MIR-72 


21:40:45.50 


58 


16:02.6 


15.68±0.09 


13.73± 0.06 


12.85± 0.04 


4.13±0.02 


3.15±0.02 


2.74±0.03 


2.88±0.03 


II 


X-ray 




MIR- 73 


21:40:45.76 


58 


15:48.9 








2.71±0.02 


7.27±0.03 


10.79±0.06 


11.29±0.07 


I 


X-ray 




MIR-74 


21:40:45.82 


58 


17:43.0 








0.36±0.00 


0.92±0.01 


1.16±0.02 


1.38±0.02 


H2 






MIR-75 


21:40:46.22 


58 


17:40.2 








0.34±0.00 


0.71±0.01 


1.20±0.02 


1.53±0.02 


H2 






MIR-76 


21:40:46.46 


58 


15:23.0 


12.81±0.03 


11.95± 0.03 


11.65± 0.02 


8.66±0.03 


5.35±0.02 


5.25±0.04 


7.74±0.07 


Trans. 


X-ray, Ha 




MIR- 77 


21:40:46.86 


58 


15:33.4 


15.3a±0.06 


13.55± 0.05 


12.63± 0.04 


5.80±0.02 


4.05±0.02 


5.16±0.04 


10.41±0.08 


Galaxy 


X-ray 




MIR- 78 


21:40:47.35 


58 


16:37.3 


17.08±. . . 


16.14±0.24 


14.63±0.10 


l.OOiO.Ol 


0.75±0.01 


0.59±0.02 


0.36±0.02 


III/S 






MIR- 79 


21:40:47.43 


58 


17:39.1 


11.53±0.03 


9.98± 0.03 


9.37± 0.02 


72.47±0.17 


46.36±0.10 


34.42±0.11 


20.21±0.07 


III/S 






MIR-8a 


21:40:48.00 


58 


15:37.6 


13.89±0.03 


12.95± 0.03 


12.67± 0.03 


4.72±0.02 


3.39±0.02 


5.77±0.04 


11.90±0.06 


I/II 


Ha, X-ray 




MIR-81 


21:40:48.52 


58 


16:37.1 


16.18±0.09 


14.95± 0.08 


14.46± 0.08 


0.71±0.01 


0.46±0.01 


0.38±0.02 


0.62±0.02 


Galaxy 


Ks 




MIR-82 


21:40:48.58 


58 


16:00.9 


15.88±0.09 


13.45± 0.04 


11.88± 0.02 


30.35±0.08 


37.01±0.08 


43.12±0.13 


61.94±0.17 167.10±0 


09 II 


JH 




MIR-83 


21:40:48.84 


58 


16:22.8 


16.02±0.12 


15.25± 0.13 


14.46± 0.08 


0.61±0.01 


0.43±0.01 


0.25±0.02 




III/S 






MIR-84 


21:40:48.87 


58 


15:10.9 


15.64±0.06 


15.04± 0.08 


14.56± 0.09 


0.77±0.01 


0.62±0.01 


0.75±0.02 


1.34±0.03 


I/II 


Ha, JKs 




MIR-85 


21:40:49.05 


58 


17:19.4 


15.79±0.08 


13.83± 0.05 


12.92± 0.02 


3.32±0.02 


2.19±0.01 


1.75±0.02 


1.18±0.02 


Trans. 


J 




MIR-86 


21:40:49.06 


58 


17:09.3 


14.14±a.03 


12.86± 0.04 


12.13± 0.02 


13.88±0.05 


14.57±0.04 


15.54±0.07 


17.83±0.07 17.99±0 


05 II 


Ha, X-ray, Ks 




MIR-87 


21:40:49.29 


58 


14:11.3 


13.52±a.a3 


12.66± 0.03 


12.35± 0.02 


3.87±0.02 


2.44±0.01 


1.66±0.02 


0.99±0.02 


III/S 






MIR-88 


21:40:49.67 


58 


17:37.3 


17.13±a.23 


14.99± 0.09 


14.23± 0.07 


1.07±0.01 


0.73±0.01 


0.60±0.02 


0.36±0.02 


III/S 






MIR-89'' 


21:40:51.07 


58 


16:24.3 


17.80±0.01 


16.70±0.01 


16.09±0.01 


O.lliO.OO 


0.08±0.00 


0.11±0.02 


0.23±0.02 


Galaxy 






MIR-90 


21:40:51.23 


58 


16:38.9 


14.23±0.03 


13.76± 0.04 


13.61± 0.04 


1.26±0.01 


0.78±0.01 


0.57±0.02 


0.53±0.02 


Trans. 


HK 




MIR-91 


21:40:51.81 


58 


16:29.1 


15.70±0.06 


14.45± 0.06 


13.93± 0.07 


1.19±0.01 


0.74±0.01 


0.60±0.02 


0.54±0.02 


Trans. 


HK 




MIR-92 


21:40:51.90 


58 


17:01.6 








0.36±0.00 


0.33±0.00 


0.20±0.02 




III/S 






MIR-93'' 


21:40:51.97 


58 


16:23.6 


19.21±0.06 


17.47±0.01 


16.39±0.01 


0.16±0.00 


O.lOiO.OO 


0.23±0.02 


0.65±0.02 


Galaxy 






MIR-94 


21:40:52.26 


58 


17:46.0 


15.68±0.05 


14.56± 0.05 


14.06± 0.06 


0.94±0.01 


0.65±0.01 


0.59±0.02 


0.48±0.03 


Trans. 






MIR-95 


21:40:52.98 


58 


18:01.1 


15.29±0.05 


14.59± 0.06 


14.22± 0.07 


0.84±0.01 


0.53±0.01 


0.54±0.03 


0.34±0.03 


III/S 






MIR-96 


21:40:53.12 


58 


14:57.4 


14.77±0.06 


14.25± 0.07 


13.95± 0.07 


0.84±0.01 


0.54±0.01 


0.38±0.02 


0.37±0.03 


Trans. 






MIR-97 


21:40:53.42 


58 


14:51.9 


15.5a±0.a6 


14.88± 0.09 


14.74± 0.11 


0.44±0.00 


0.31±0.01 


0.23±0.01 


0.20±0.02 


Trans. 






MIR-98 


21:40:53.68 


58 


16:41.2 








0.25±0.00 


0.71±0.01 


0.79±0.02 


0.83±0.02 


H2 







Table 1 — Continued 



Source 
SF038 



<52000 



^3.6 



mjy 



F4.: 



mJy 



mJy 



Fs.o 
mJy 



F24 
mJy 



-99 


21:40:55.68 


58:15:50.8 








0.36±0.00 


0.14±0.00 


2.69±0.03 


9.65±0.05 


-100 


21:40:55.87 


58:14:16.3 


14.18±0.03 


13.61± 0.04 


13.39± 0.04 




0.87±0.02 




0.45±0.03 


-101 


21:40:55.88 


58:15:34.8 








0.25±0.00 


0.32±0.00 


0.41±0.02 


0.71±0.02 


-102 


21:40:55.93 


58:17:07.8 


14.45±0.03 


13.19± 0.04 


12.76± 0.03 


3.21±0.02 


2.01±0.01 


1.39±0.02 


0.99±0.02 


-103 


21:40:57.04 


58:16:56.8 


16.61±0.13 


15.92± 0.21 


16.22± 0.00 


0.42±0.01 


0.34±0.00 


0.40±0.02 


1.19±0.03 


-104 


21:40:57.06 


58:16:29.0 


14.73±0.03 


13.00± 0.03 


12.24± 0.02 


6.10±0.02 


4.11±0.02 


2.87±0.03 


1.96±0.03 


-105 


21:40:57.33 


58:16:41.7 


16.64±. . . 


15.41±. . . 


14.55± 0.13 


0.79±0.01 


1.18±0.01 


2.90±0.03 


6.34±0.04 


-106 


21:40:57.50 


58:14:43.6 


15.00±0.04 


14.66± 0.07 


14.34± 0.09 




0.43±0.01 


0.25±0.03 


0.15±0.02 


-107 


21:40:58.03 


58:15:12.1 


15.06±0.05 


14.50± 0.07 


14.12± 0.08 


0.70±0.01 


0.46±0.01 


0.30±0.02 


0.21±0.02 


-108 


21:40:58.45 


58:16:20.5 


13.11±0.03 


12.44± 0.03 


12.02± 0.03 


5.93±0.03 


3.77±0.02 


3.13±0.03 


2.52±0.03 


-109 


21:40:58.53 


58:15:00.3 


12.15±0.03 


11.26± 0.03 


10.96± 0.02 




8.34±0.03 


5.73±0.05 


3.46±0.03 


-110 


21:40:59.71 


58:17:33.3 


12.69±0.02 


11.25± 0.03 


10.59± 0.02 


24.86±0.17 




10.77±0.13 





I/II 

III/S 
Galaxy 

Trans. 
Galaxy 

Trans. 

H2 

III/S 
Trans. 
Trans. 

III/S 
III/S 



'^JHKs photometry from lNisini et al.l )2001l ) 
^JHKs photometry from lBeltran et al.l 1200911 
■^Variable in J/H/K^ 



to 

00 
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Table 2. Coordinates and flux densities of NIR & IRAC 3.6/4.5 fim PMS sources in SFO 35 



Source 


"2000 


^2000 


J 


H 


K, 


F3.6 


F4.5 


SF038 












mjy 


mJy 


NIR-1 


21:40:34.36 


58:16:19.2 


20.88±. . . 


18.12±0.03 


16.68±0.02 


15.73±0.03 


15.38±0.03 


NIR-2 


21:40:36.46 


58:16:28.7 


20.84±. . . 


19.31±0.10 


17.14±0.03 


15.23±0.02 


14.75±0.02 


NIR-3 


21:40:37.57 


58:16:01.2 


20.80±. . . 


19.98±. .. 


17.35±... 


15.22±0.03 


14.29±0.02 


NIR-4 


21:40:38.17 


58:15:32.7 


20.87±. . . 


18.32±0.03 


16.72±0.02 


14.84±0.02 


14.59±0.02 


NIR-5 


21:40:39.53 


58:15:46.9 








14.22±0.01 


13.72±0.01 


NIR-6 


21:40:39.65 


58:15:06.6 


16.59±0.01 


15.88±0.01 


15.52±0.01 


14.07±0.01 


13.82±0.01 


NIR- 7 


21:40:42.81 


58:15:45.3 


16.29±0.00 


15.38±0.00 


14.77±0.00 


13.91±0.01 


13.27±0.01 


NIR-8 


21:40:44.15 


58:15:49.6 








15.78±0.04 


13.83±0.01 


NIR-9 


21:40:44.57 


58:15:51.6 


15.97±0.00 


15.23±0.00 


14.71±0.00 


14.28±0.02 


13.68±0.01 


NIR- 10 


21:40:45.85 


58:15:40.6 


18.70±0.03 


17.09±0.01 


16.01±0.01 


14.68±0.02 


14.18±0.02 


NIR- 11 


21:40:48.28 


58:16:18.8 


18.17±0.02 


17.16±0.01 


16.40±0.02 


15.63±0.03 


15.24±0.03 


NIR- 12 


21:40:50.55 


58:16:01.8 


20.84±. . . 


19.65±0.10 


17.55±0.03 


15.65±0.03 


15.29±0.03 


NIR- 13 


21:40:51.14 


58:17:01.5 


20.86±. . . 


18.74±0.05 


17.23±0.04 


15.99±0.03 


15.68±0.04 
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Table 3. BVRI photometry of MIR sources in SFO 38 



Source 


B 


V 


R 


/ 


Variability 


MIR-1 


19.82±0.010 


18.24±0.005 


17.18±0.004 


16.00±0.003 




MIR-2 


20.48±0.015 


18.88±0.006 


17.88±0.008 


16.84±0.005 




MIR-3 


21.58±0.021 


19.58±0.007 


18.37±0.006 


17.14±0.004 




MIR-4 


O.OOiO.OOO 


20.38±0.016 


19.15±0.007 


17.73±0.006 




MIR-5 


21.11±0.015 


19.32±0.007 


17.77±0.004 


16.37±0.002 


BVI 


MIR-7 


21.59±0.021 


19.39±0.006 


17.90±0.004 


16.38±0.002 




MIR-8 


17.74±0.001 


15.74±0.001 


14.55±0.001 


13.35±0.001 




MIR-9 


19.95±0.017 


18.48±0.004 


17.60±0.005 


16.62±0.004 




MIR-10 


17.28±0.001 


16.07±0.001 


15.29±0.001 


14.49±0.001 




MIR-11 


O.OOiO.OOO 


19.52±0.011 


17.94±0.007 


16.19±0.005 




MIR- 14 


17.67±0.001 


16.24±0.001 


15.38±0.001 


14.53±0.001 




MIR-15 


16.89±0.001 


15.26±0.001 


14.33±0.001 


13.44±0.001 




MIR-18 


19.62±0.006 


17.58±0.003 


16.30±0.002 


14.98±0.009 




MIR-26 


17.34±0.001 


16.03±0.001 


15.20±0.001 


14.30±0.001 




MIR-29 


20.58±0.011 


18.43±0.005 


16.94±0.002 


15.42±0.002 




MIR-31 


17.43±0.002 


15.71±0.001 


14.60±0.001 


13.52±0.001 




MIR-32 


21.15±0.015 


19.29±0.006 


18.01±0.007 


16.25±0.002 


BVI 


MIR-43 


19.41±0.004 


17.49±0.002 


16.21±0.001 


14.73±0.001 


BVI 


MIR-44 


19.36±0.006 


17.65±0.003 


16.57±0.001 


15.51±0.002 




MIR-45 


21.41±0.018 


19.45±0.007 


18.10±0.012 


16.33±0.002 


BVI 


MIR-47 


O.OOiO.OOO 


19.86±0.009 


18.76±0.005 


17.61±0.005 




MIR-49 


20.91±0.012 


18.81±0.004 


17.04±0.003 


15.69±0.001 


BVI 


MIR-52 


19.54±0.008 


18.32±0.005 


17.52±0.003 


16.64±0.005 




MIR-58 


20.13±0.007 


18.75±0.003 


17.86±0.005 


16.91±0.003 


B 


MIR-60 


20.79±0.011 


19.13±0.005 


18.07±0.005 


17.06±0.003 




MIR-64 


O.OOiO.OOO 


20.49±0.014 


18.77±0.016 


17.09±0.006 




MIR-76 


19.08±0.006 


17.16±0.002 


15.90±0.002 


14.55±0.001 




MIR-79 


21.41±0.024 


18.56±0.004 


16.52±0.002 


14.55±0.001 




MIR-80 


21.69±0.023 


19.81±0.008 


17.90±0.006 


16.00±0.004 


V 


MIR-87 


19.72±0.005 


17.72±0.002 


16.50±0.002 


15.28±0.002 




MIR-90 


18.15±0.002 


16.85±0.002 


16.04±0.001 


15.26±0.002 




MIR-95 


O.OOiO.OOO 


20.16±0.012 


18.71±0.009 


17.34±0.004 




MIR-96 


19.23±0.008 


17.85±0.003 


16.94±0.003 


16.02±0.002 




MIR-97 


20.06±0.007 


18.63±0.005 


17.76±0.005 


16.83±0.004 




MIR-100 


18.31±0.003 


16.86±0.002 


16.03±0.001 


15.28±0.002 




MIR-106 


18.87±0.004 


17.75±0.003 


16.96±0.003 


16.21±0.002 




MIR-107 


19.50±0.004 


18.12±0.002 


17.22±0.003 


16.27±0.002 




MIR-108 


18.94±0.005 


17.31±0.004 


16.13±0.003 


14.86±0.001 




MIR-109 


18.17±0.003 


16.17±0.001 


15.00±0.001 


13.86±0.001 




MIR-110 


O.OOiO.OOO 


19.78±0.010 


17.76±0.004 


15.76±0.003 


V 


BRC38 1 


O.OOiO.OOO 


20.05±0.017 


18.24±0.009 


16.45±0.003 




BRC38 17 


O.OOiO.OOO 


20.51±0.017 


18.65±0.005 


16.92±0.006 


V 


GFG 64 


19.01±0.007 


17.48±0.003 


16.19±0.001 


15.18±0.001 




IPHAS 


O.OOiO.OOO 


20.28±0.011 


18.28±0.009 


16.50±0.004 
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Table 4. Spectral Classification of YSOs 



ID 


Alt. ID 


Date 


(A) 


Sp. Type 
(isubtype) 


(K) 


Type 


Ho[10%] 
(km s-i) 


log Mac 

(M0yr-i) 


aiRAC 


E(V-I) 


Ay 


E(V-I)/E(B-V) 


BRC 38 1 




5/10/08 
10/10/09 


-12.93 
-16.24 


M3±l 


3470 


CTTS 


312.23 
205.29 


-9.86 
-10.90 




1.13 


2.81 




MIR-5 


BRC 38 2 


5/10/08 
10/10/09 


-85.23 
-106.1 


Ml±l 


3720 


CTTS 


359.79 
244.53 


-9.40 
-10.52 




0.99 


2.47 


3.19 


MIR-7 


BRC 38 3 


6/10/08 
10/10/09 


-19.0 

-18.75 


Ml±2 


3720 


CTTS 


476.25 
463.50 


-8.27 
-8.39 


-0.87 


1.05 


2.61 


1.46 


MIR- 11 


BRC 38 4 


10/09/08 
18/07/09 


-7.30 
-11.29 


Ml±2 


3720 


WTTS 






0.02 


1.37 


3.41 




MIR-29 


BRC 38 5 


6/09/08 
3/11/08 
24/08/09 
9/11/09 


-3.42 
-8.84 
-1.89 
-2.74 


Ml±2 


3720 


WTTS 


469.90 


-8.33 


-1.74 


1.05 


2.61 


1.57 


MIR-31 


BRC 38 6 


6/09/08 

24/08/09 

9/10/09 


-53.47 
-57.32 
-60.5 


K2±l 


4900 


CTTS 


384.43 
442.02 
382.26 


-9.16 
-8.60 
-9.18 


-1.60 


1.18 


2.94 


1.42 


MIR-32 


BRC 38 7 


10/09/08 
12/09/09 


-33.12 

-15.78 


M2±l 


3580 


CTTS 


368.03 
220.59 


-9.32 
-10.75 


-1.03 


0.90 


2.24 


2.65 


MIR-43 


CFG 62 


18/07/09 


-3.516 


Ml±l 


3720 


WTTS 






-2.61 


0.80 


1.99 


1.82 


MIR-45 


BRC 38 9 


6/09/08 
3/11/08 
24/08/09 


-37.43 
-41.22 
-56.33 


K4±l 


4590 


CTTS 


281.09 
389.74 
273.39 


-10.16 

-9.11 

-10.24 


-0.18 


1.97 


4.91 


2.21 


MIR-49 


BRC 38 10 


4/11/08 
9/10/09 


-13.5 
-30 


M0±1 


3850 


CTTS 


430.70 
450.55 


-8.71 
-8.52 


-1.37 


1.32 


3.29 


1.91 


MIR-64 


BRC 38 11 


10/09/08 
12/09/09 


-31.46 
-25.06 


K5±l 


4350 


CTTS 


239.20 


-10.57 


0.10 


2.04 


5.08 




MIR-76 


CFG 81 


18/07/09 
10/10/09 


-3.261 
-3.135 


K7±l 


4060 


WTTS 






-1.07 


1.01 


2.51 


1.87 


MIR-80 


BRC 38 12 


24/08/09 


-8.404 


M4±l 


3370 


WTTS 






0.31 


0.95 


2.37 


3.39 


BRC 38 17 




11/11/09 


-119. 


M3±2 


3470 


CTTS 


205.29 


-10.90 




1.12 


2.79 





CO 
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Table 5. Parameters derived from SED modeling using axisymmetric radiation transfer models 
for the candidate YSOs. A^ refers to the foreground extinction towards the source. 



Source 


n 


M, 


Mdisk 




Mcnv 




-'"cnv 


Mdisk 


L 


A. 


Age 




K 


Mo 


Mo 




Mo 




Mo yr-i 


Mo yr-i 


Lo 




Myr 


MIR- 11 


4442 


1.26 


6.44x10- 


-3 


6.39x10- 


-5 


0.0 


2.10x10-9 


1.27 


5.3 


3.2 


MIR-29 


4837 


1.55 


1.66x10- 


-3 


6.13x10- 


-7 


0.0 


1.16x10-9 


1.67 


4.4 


5.8 


MIR-31 


5923 


1.88 


1.33x10- 


-4 


7.04x10- 


-8 


0.0 


2.41x10-9 


8.92 


4.2 


6.7 


MIR-34 


4523 


1.53 


7.64x10- 


-5 


2.56x10- 


-3 


1.21x10-^ 


1.28x10-10 


5.8 


3.2 


0.69 


MIR-36 


5413 


3.29 


4.37x10- 


-3 


1.85x10- 


-8 


0.0 


4.80x10-*^ 


29.1 


49.6 


1.3 


MIR-38 


7854 


1.98 


3.10x10- 


-5 


2.99x10- 


-6 


0.0 


5.49x10-" 


11.9 


23.0 


8.2 


MIR-45 


9960 


2.32 


4.23x10- 


-5 


1.70x10- 


-8 


0.0 


1.37x10-9 


27.0 


8.3 


9.4 


MIR-48 


7554 


5.14 


1.71x10- 


-3 


21.4x10- 


-2 


2.23x10-^ 


2.38x10-9 


280 


45.9 


0.49 


MIR-50 


4643 


5.97 


7.67x10- 


-2 


40.8 




3.82x10-4 


2.08x10-^ 


196.7 


44.4 


0.1 


MIR-54 


4290 


1.50 


2.46x10- 


-2 


9.3 




7.13x10-'^ 


2.98x10-6 


33.4 


19.6 


0.17 


MIR-55 


4989 


3.07 


7.86x10- 


-4 


1.44x10- 


-3 


1.97x10-9 


1.14x10-9 


14.4 


52.3 


0.91 


MIR-59 


5218 


4.27 


9.99x10- 


-3 


7.37x10- 


-3 


2.95x10-^ 


2.00x10-^ 


60.5 


87.4 


0.53 


MIR-82 


11320 


2.77 


5.85x10- 


-5 


1.18x10- 


-6 


0.0 


1.40x10-19 


57.2 


21.0 


4.8 


MIR-86 


4594 


1.55 


8.18x10- 


-4 


3.63x10- 


-7 


0.0 


4.50x10-10 


2.59 


9.7 


1.8 



